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ABSTRACT 


A bubble transport equation was developed: and. a: formai_ solution 
obtained for the bubble distribution finction.. The: distribution- 
function obtained depends upon. the model used. for: bubble-drag;.gas- 
diffusion, and acean circulation.. As a. specific: application, .solutions: 
were obtained for a one—dimensional,. quiescemnt,.steady state-ocean 
with a bottom source at I97 metenrs.. The. characteristié: equations - 
were integrated numerically to test the accuracy of* approximate. 
analytical solutions. The calculated bubhle densities were’ compared 
with experimental data available for this problem.. When gas diffusion 
is neglected, the theoretical and: experimental distributions are.in- 


Pemeral agreement. 
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Ilo JINTRORUELION 


Studies have been made by many researchers regarding the mechanisms 
of gas bubble dynamics in a flutd from both the theoretical. [Ref:.1]_ 
and the experimental [Refs. 2, 3, and 4] aspects.. The mechanisms: of_ 
gas bubble formation, transport, - extinction. are important. as. gas 
bubbles can effect the propagation of sound in. the ocean by the under-- 
water scattering of sound beams.. Theoretical studies have explained 
how a bubble varies in size as it rises in. a quiescent fluid:. The 
phenomena involved includes shrinkage from gas diffusion. versus 
expansion from decreasing hydrostatic pressure and shrinkage from 
surface tension for very small bubbles.. However,. there: does not appear 
to be at present any theory explaining how the distribution of bubbles 
generated by a bottom or distributed source in. the ocean.cam be: expected 
to vary with depth or initial bubble radius.. 

This paper outlines the theory for determining theoretically the 
distributions of bubbles generated by a bottom source in the ocean 
[Ref. 5]. First, a linear bubble transport equation is developed and 
then it is solved by the method of characteristics... A model fora one-- 
dimensional steady state ocean is specified for application of this 
theory. Depending on the assumptions made regarding the dynamics 
involved in the bubble motion, such as the importance of gas diffusion 
or surface tension, various semi-analytical solutions can.be obtained 
for the distribution of bubbles a a. bottom source.. These: distri-- 
butions are then compared to available experimental data.. In particular,. 
the distributions from a bottom source located at 197 meters. were calcu-- 
lated using various approximations and the results were compared with 


the experimental data in Ref. 4.. 
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i. WE BUBBLE TRANSPORT EOUATION 


The transport equation is basically a balance equation [LRefs.. 6 
and 7]. It was desired to develop a transport equation. for the. 
description of the distribution of bubbles as @ function of position,. 
velocity, size (bubble radius) and time in a circulating field’ that 
may contain sources or sinks. A distribution. function that describes 
the bubble population can be defined and if the bubbles satisfy the 
following general dssumptions, this distYrbution will satisfy a: 
Boltzmann-type transport equation: 

(1) The bubbles are completely independent of one another;. 

(2) If bubble interaction occurs, no more than. two bubbles will 
interact at any one time. 

Since bubbles in the ocean exact no force on one another andthe 
history of one bubble does not affect the history of another bubble, 
the two above assumptions hold and a transport. equation.was found.. 

et w(r,v,2,t }d3rd3vag is the number of bubbles at time t: with 
radius dg about 2, velocity d3y about Vv, and position d3r about Yr, 
then at time t' = t + dt, iGo ag tt) d3r"d3v'dg' is the number of 
Pebbles at time t + dt with radius dg' bout 2',. velocity d3v' about wv’, 


. 'e -> — ->,. : 
and position d3r' about r', where r', v', and g" are: 


r' =P + vdt 
vw’ = ¥ + Adt (1) 


Cl = fete ee 


The above Taylor series expansions about time t neglects all terms. of 


order two and higher. The changes in position, velocity, and radius 








are the result of nonzero velocity ¥ = dr/dt,. acceleration a =-dv/dt,. 
and radius change rate y = dg/dt, respectively.. 

Since the transport equation is a balance equation the- number. of 
bubbles in a volume element d/t" = @r'atvi ae" at- time: t- + dtyz.must 
equal the number of bubbles that were in. the volume d’t =: d?rd?vag 
at time t plus or minus the bubbles that were generated or removed by 
any internal sources or sinks. Thus an equation. for the. conservation~ 


of bubbles is written as: 


yr" 2v" oh stedt)}d7c° = W(D, v2 »t)d7r = S(rv 2 ee Tt (2) 


where S(r,v,2,t}da7r is the rate at which bubbles were introduced into 
the volume element d’t about (T,752) by distributed® sources... 


: i a J F ee >-->-- 
wr? vv" 2" ttt) in. a Taylor series about. w. ==(r,v,2,t) 


Expanding w' 


and keeping only terms through the first order yields:: 


vi = wv t+ Cav/fat)dt + (7 + V_w)dt + (a- V_widt: + v(dyp/ot)dt- (3) 
where 
v (>) = 942)/dx4.01 (7d a --)/ 92 
vata} = at) /ov_ > Sitae acwenene em =: 


The volume element d’rt" was transformed into terms: of d’t by using 


a Jacobian. The results of this transformation yielded. 


g7qt = MEV! 2") ge (4) 
>> 
a(r,v,2) 
with 
> > > 
MO ae Te (y~ a)dt + (av/aeyat.. Cs) 
s(P,V,2) | 





Substituting (3), (4), and (5) into (2),. dividing by d/’tdt. 


and taking the limit dt > 0 ytelds the bubble transport equation 


(eRe fee oul: 
wi >. 2 a tie Ye. > ie | 
ae an woe ree eu at a ce as x0 yeas Oe 


It should be noted that this transport equation is: different: from 
7 ss 
the usual transport equation encountered because the term (V sa + dv/dt) 
e - - oe ve ana . 
is generally nonzero as the bubble radius change rate,. v(r,v3;2,t),. and 
® er ae 2 - = e@. - - - 

acceleration, a(r,v,2£,t), are functions of bubble radius and- bubble 
velocity respectively. 

The distributed sources or sinks are contained in the term 


>> 
S(r,v,2£,t). This term can also include the rate at which. bubbles 


S 
are introduced into d’t and removed from d’t by scattering off small 
scale turbulent eddies or by bubble-bubble interactions.. It was. assumed. 
the S was independent of or proportional to {.. This simplication 


yielded a first order, partial differential bubble transport equation.. 





ETL. FORMAL SOLUTEON TO THE. TRANSPORT EOURie) 


The bubble transport Eq. (6) in the form of a first order: partial 
differential equation was solved using the method of characteristics 
[Ref. 5]. This was equivalent ta replacing the transport. equation 
with a set of simultaneous, first order,. ordinary differential equations.. 

The characteristic equations were determined to. be the bubbie 


dynamic equations: 


dp/dt = 4 mit =r = (x.,7 2) (7a) 
Go oO Oo © ©: 

> _> >>. ie 

dg/dt = y A cra) a) Cia} 
oO oO 

dy/dt = S - Ly wt v= mo (7d) 


The model chosen determined the acceleration a(r,v32,t) and rate of. 
bubble radius change: v\(r,v,2,t),-thus establishing the family of. char-- 
acteristic curves. This family was parameterized by the set of-initial 

° e > cated > ae eT eet wmen: 2 ° ° ° 
Bemcitions (ri,v.5k_)- andy (2 ,Vv_,2.) was sthe initial distribution=. 

roo OO =e 

The general form of the distribution function was obtained by 


integrating Eq. (7d) and making use again of a Jacobian.. The formal 


Borution for the bubble distribution is’ fRef.. 5]: 
Fc 


ty 


w(r,v,2,t)d7r = pCR ev 2. ot) SUE s(piyy oe" ph jda7 matt 
te 

The integral here is a line integral along a characteristic: curve as 

are all the integrals in this thesis. Equation (8) is equivalent to 

Saying that the number of bubbles in. d7~ about (457,52) at. time: tis. 

the number of bubbles that were originally in dae. plus or mihus those 


added or subtracted by sources or sinks while proceeding along the 


characteristic curve. 





IV. THE MODEL FOR THE SOLUTION OF THE BUBBLE TRANSPORT EQUATION 


- Before solutions to Eq. (6) could be obtained explicitly ,. a: model. 

of the upper ocean had to be specified.. This model must provide: explicit 
e _ eee ae : 

expressions for the acceleration a(r,v,2,t),. the rate of change of. 


=e 


bubble radius v(%,7,2,t); the bubble distribution S(v,¥752,t) and_the 
-> 
circulation field V(v,t). A Cartesian coordinate. system was: chosen 


with the origin located at the surface of the ocean. and. the z-axis 


vertical and positive upward. Thus all depths were taken to be. negative.. 


A. THE BUBBLE RADIUS CHANGE RATE 

The rate of change of bubble radius for bubbles circulating im the 
ocean is caused by either gas diffusion, @ change in the hydrostatic. 
pressure, surface tension, or a combination of these phenomena.. This 
interaction of effects was analyzed by assuming the bubbles to be. 
spherical containing an ideal gas. From the ideal gas law 

ee Se (9) 

where R is the gas constant, T is the temperature, P is the pressure, 
and V is the bubble volume. By taking the time derivative of Eq. (9) 


and assuming isothermal changes, an expression for yv,. the change. in 


bubble radius was obtained: 


Gi Se Gi, Ge 
a ae 0 
v= at Pier sz sy 
The term - = contains the gas diffusion effects while the effects of 


the change in hydrostatic pressure and surface tension are: contained 


: 1 dP 
in the term 5 ap z 





The hydrostatic pressure P is the sum of the atmospheric’ pressure,. 
cai and the fluid pressure, “P82 + The surface tension pressure is 
2y/2% , where  , is the density of sea water, = 1.03 x 10° Kg/meter®, . 


and y is the surface tension of sea water,. approximately: equal.to.0.74- 


newtons/meter. Therefore 


= — g 
P Se pez + 27 (11) 
and 
dP ay ak 
| Ja 12 
at Poe ~ 22 at er 


By using a simple model for gas diffusion [Ref.. 8] the: following. 


expression was obtained for dn/dt 


a — éturg (Pp - EB? (13) 


where § is a constant. A value for § had been. determined from experi-- 

mental data to be about 1076 meters/second... It was assumed for. this 

model that the dissolved gas concentration gradient was. confined to 

a thin shell at the air-water boundary because of the bubble's speed. 
An expression for y was obtained by substituting Eqs..(11),.(12), 


mat ls) into Eq. (90) [Reft. 54. 


ve A = [v 22 + B6RT(ze — r)I/[3(D-2)e + 2P] (14) 
where 
T = 2y/p = 1.47 x 10 > meters? 
and 
» =P 46 ¢ = 10 meters. 


oro 


If 2 is greater than 30 microns (I micron (yu) =10°* meters), the 


surface tension terms can be ignored,. i.e., T = 0. 


10 





B. THE BUBBLE EQUATION OF MOTION AND ACCELERATION 

A bubble that is in motion with respect to a stationary reference 
system will experience several forces acting on it.. First. there:is- 
the force of the water acting on the bubble to pull it-along:.. Ifthe 
bubble were completely entrained by the fluid so that: the bubbie’s: 
velocity (7) were equal to the circulation Field velocity (V),.the 
entrainment force would be o adV/dt. where g is the: volume of: the: 
bubble. Complete entrainment can. occur only: if @ is’ equal to-p... 
Since this is rarely the case, the bubble will maintain a relative- 
velocity u = y — V with respect to the water.. This relative velocity 
produces a drag force, ee acting on the bubble to slow it. down.. 

As fie Wibpic meves throtgh the water,. the bubble: pulls water. 
along with it. The amount of water pulled along with it is equal_to 
some fraction, 8, of the bubble's volume.. From experimental data-8 is 
usually taken as 0.9. This phenomena is: called the mass adhering force 
and is equal to -(Bop du/at).. The bubble will also: experience: a:net. 
bouyant force which is equal to glo | —o) Z,. where g is the gravita- 
tional acceleration. 

By using Newton's second law these forces were: combined to: yield 


fe 
an expression for the acceleration a L[Ref.. 5]: 


= 
FE i a 
a i ile gs ae OS ge Crs) 
op (Bte/e.) (B+o/p ) dt (B+o/o .) 
where 

ian oO ee eG 

ee -V ~ 16 

One ot ihe ay (16) 


—_-- 
In order to obtain an explicit expression for the: acceleration: a,. 


>- 
it was necessary to derive an expression for the drag force Pot Various 


JE 





models for the drag force were considered, with the models being 
classified by using the Reynold"s number R, = 2hup /n where n is- 
the viscosity of water (n = eae Newton-second/meter?).. 

For bubbles with a radius less than or equal to 100 microns. 
(2 < 100 ), the Reynold's number fs less thar one (R, <1) and-the 
drag is viscous according to Stokes Law: 


pe 


ED Sat 


- Ky Zu (8 
where K = 67 in the presence of surface active material and: 4" in- the 
absence of these surface active materials.. 


The drag was considered to be both viscous: and. turbulent- if: 


moon <£< 1500p0r 1 < Ro < 700. Thus the drag force was 


> > ae 
FD : | Geakes i Erurb sa 
ae 3 . * . eae a 
Here Fo, peg Was given as in Eq. (17) with 60s K s 10% sand Frurb? Which 


epise from the energy dissipated in the turbutence of the wake>.hadethe 
form: 


—-~ 


E = comeing ie (19) 
Turb ae fore 


where Ke is a drag coefficient and Sy is the hubble surface. area’ 
covered by the turbulent wake (s, = 72%). If there were no surface 


ie at e.- 
active materials present PS would equal F > Waa K == 12124 


Stokes. 
Since the oceans contain surface active materials,. the drag- force 
contains both a viscous drag and turbulent drag term. By noting that 


p/P | << 1 and thata = 4237/3 ror Dub les Wit Vi oOOl Wy Gi ec Colma. 


arion Can be written as@ Eker. Si 


B= mauler oa re ee ee ee (20 
2 Q B B 
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a = 2 aie, = 8.75 x 10 © meters/sec? 
2 Bo 
O 
and 
ake Sy : 
Sa = 4S a LO 
ume? 


If the bubbles were generated by a surface source: such as rain, 
wave action, etc., the bubble radii are generally less than 1507 y 
[Refs. 3 and 9] and the drag force is predominately-viscous.. In-this- 
case —€ should be set equal to Zera.. 

If the bubbles originate at a bottom source,. the radii are inthe 
range of 500 yto 700yu [Ref. 4] and the drag should be considered to be 


predominately turbulent with a set equal to zero.. 


tee BUBBLE SOURCES 

The bubbles in the ocean might physically be considered to have 
originated from three types of sources.. The first type included 
bubbles that originated near the surface and. were caused by wave- 
action, rain or snow. The distribution of bubbles originating-from 
these phenomena was analyzed in Ref. 10.. 

The second type of bubbles are these nucleated by distributed 
particles of material as well as bubbles that originated from biological 
activity such as photosythesis or bacteria and are distributed through- 
out the medium. 

The third type of bubbles includes those bubbles generated on the 
floor of the ocean. These bubbles are often sediment-—initiated bubbles 
and measurements indicate these bubbles to he between 450 p-- 800 yin 
radius with a peak around 550y[Ref. 4J. It was this category of 


bubbles for which the distribution function was evaluated. It was 


1S 





considered to be mathematically convenient to introduce the bubbles 
via boundary conditions rather than to inelude them in S(r,v,2,0).. 
y Beem represented the bubble distribution originating ir 
the bottom sediment. 

The model for the circulation field was chosen to be @ quiescent 
fluid dictating thet V=0. For this model the bubbles could: rise. 


only via buoyancy and the characteristic curves are vertical upward.. 


14 





Va SheESULiS 


For this thesis the bubble transport theory was applied: to a: 
one-dimensional quiescent, steady state ocean with a bottom source 
at 197 meters. The Seton and auprorinete coluion=] eee 
The characteristic equations were integrated numerically to test 
the accuracy of the approximate analytical solutions.. The- calculated 
bubble distributions were compared with. experimental. data. 

In thas steady state model the velocities are vertical upward,. 
and all functions depend spatically on. depth only.. The: bubble transport- 


equation (Eq. 6) can be written as: 


CO ie eee: = 
Weae 12 dv TUG zo oe (21) 


Here $¢(z,v,2) is the distribution function,. S is the. distributed source. 
munetLOn, v is the vertical velotity,. v is the rate of change: of- bubble 
radius, and a = ga/av + av/ar.. 

If z is chosen as the independent parameter,. the characteristic 


equations are [Ref. 5]: 


Coe i ao g 

aie ge Q Nr BV see) 

dQ ov _ 22 + S&RT( ze - T)/v 

eee 3(D-z)2 + 27 See) 
and 

ag 1 fe | 

@z~ 7 =v tH? tEFe pez) 


Since the characteristic curves are vertically upward‘ along~ the. z=axis,. 
the characteristics or boundary conditions are v = Vik = 2 ,.and 


o = OZ sv 2k) at 7 = Z 0: 


iS 





Equation (22c} was integrated ta yield 


Z 
2 = (ZV oko) o[- | raw" 2" azt/v"| 
a Zee) oO 


+ exp [ - £2" 0" 2" az" /ve | SUZ Nk id ze a Gay) 


tT 
Zi 


O 
with the integration taking place along the characteristic curves. 
{z,vCz32 4V, 2) 20252 5 oh )5 Zz < 0}. )Bym@isang the Jacobian Equal 2s) 
was rewritten as [Ref. 5]: 
vo(z,v,2)dvdg = Vigl2 sv 2 )dv de. 
Z 
+ Sh7' VE az ay dc". (24) 


ZO 
Equation (24) is interpreted as @ conservation equation which states 


that the flux of bubbles, vo, at the point (z,v,2) is equal. to. the 
flux of bubbles, Vio¢o? at ae point (2 ,V .2)) plus (orn minus) any 
bubbles added (or substracted) by distributed sources along the 
characteristic curve. 

Prior to the calculation of the distribution function in either 
Eq. (23) or Eq. (24) it was necessary to establish a mapping between. 
meants (z.,v.,2) and t(z",v' .%°') on the characteristic curve.. This was 
done by integrating simultaneously the characteristic equations (Eq. 
a) touobtain vizyes 2" wid) ) ei 0(z pee 9!) for weieious. 
Z",v', and g' and substituting these relationships into Eq. (23). 
Since a/2g2 >> 1, it was very time consuming to integrate. Eq.. (22a) 
numerically. Referring to Fig. i it is seen that an. appropriate 
approximation is to equate the velocity and the terminal velocity.. 


This permits semi-analytical solutions.. 
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This approximation can be verified by rewriting Eq. (22a) as 


ipek. 5] 
= = > Co View) (ve + v") (25) 
where Yo 
-l+ [1 + 4(v_./v_,)71 
- Tlie Te 
ae ee et (26) 
2( V/V”) 
d 
os EACLE Mpg /%p9)2T (27) 
. | 
iy 
ee ay 


is the terminal velocity in the case where viscous drag dominates 


(g < 100 y). 


[z= 


(28b) 


i /2 


g 
i fe 


is the terminal welocity in the case where turbulent drag dominates. 


wm 


(g > 400 y)- Thus 


Vv. for ge< 100 


TL —. 
Yn = a (29) 
> 
Vo Lor Q > 400 


Equation (25) was integrated using fourth rank optional. Runge Kutta 
[Ref. 11]. These numerical results are presented graphically in Fig. l. 
It is observed that when a bubble is released from a depth of 20.00 
meters with zero velocity, It will reach terminal velocity after. having 
traveled a distance of only several bubble radii... Thus the: approximation 
‘that v = Ver is considered to be a valid approximation for this model. 
This approximation is convenient as Eq. (22a) is very difficult to 


integrate numerically. 


ee 








By uSing the approximation v = Vie the bubble distribution function 


can be written approximately as: 


$(z,v,L} = d(v - Vir) C252) 


where 


co 


(22) r= | GZ ved dy 


—_O 


(30) 


(31) 


ms the "depth-radius" distribution fumciion. By substitutinmggeq. ~Ga 


into Eq. (24) and integrating over v, an expression for ¢ was: determined -_ 





Bemiref. 5]: 
Vv AR 
_ ale O 
(2,2) = s— oCz 220) — 
al f, 
where sa 
Si 0!) 2  Sén" " aa" 
and = 
ae 
=j—/ d ee 
de. a 2 | eve 


(22) 


(33a) 


(33) 


In order to calculate 6(z.¢), it was necessary to determine 92'/92 


and £'(z';z,&%) either analytically or numerically. Note that 


AL*'/dX = 1/(dR/dL"'), and that the approximate differential equation 


that dYd2' satisfies is obtained by setting v = Ver 


taking the partial derivative with respect to 2'[Ref.. 5]: 


as °((3) [-(CE-]] /bo 3] 
[2 (22 )fe+-SG)/a [If / fee 4] 








where 
= 
a 
ee aa b not ae Ve 
yon vo ee 
; “e 
Furthermore, with the approximation v = v,, then = Rens and Eq. 


ars 


in Eq..(22b) au@ 


(34a) 


(34b ) 


(22b) 





becomes 


dy g2 + 36RT(z2 — [)/v 
= I (35) 
dz” 3(D-z 2 + 2F 


Equations (34) and (35) can be integrated numerically using optimal. 
fourth rank Runge Kutta [Ref. IL]. These results car be substituted 


meee bq. (32) to obtaim®= ¢. for varicus sources.. Since vy =v o15) 


a Tg 
a good approximation, it is expected that the approximation &==%, will. 
also be very good. 

To obtain approximate analytical results for comparisor.to- the 
more exact numerical results outlined above,. the approximations: of: 
Eq. (29) can be used. Neglecting surface tension,.l = 0 is aiso a. 


Valid approximation for % »30 y. Them Eq.. (35) can be integrated 


emalytically [Ref. 5}. For wv = Very alee = 0 


73: 
D-z , «4 3 adRT Zz? oi2 
Ca oe o7*® ZF} = [—— QL! — a : 
7b 232 32* } E + 





2 (s/b) D-z (36a) 
Te 


D-z'" , ORT {p- ' 
Peer (2,2 5¢') = ( 3 A! a fa) i ! 
= _ ((g/8)6) 1/2 \ 


ay l/2 Cove 
- 3D € - (F ; )| a (37a) 


ime inverse transformation are 


mor vy = Vv and, [= © 








NEN 





oo ea =a.) (e6b) 
and 


oe. ee 37b 
respectively. Then 


og? D-z g2 
Se ame ey 


TS 





and 








1/2 1/2 
ok* ey baz : 
-— = (2-2, [2 | (37¢) 


respectively. If gas diffusion can. be neglected,. 6 = 0,then Eqs..(36) 


and (37) reduce to: 


, r7 3 
ee ~z" ; 
8 = (38a) 


; = e 


Figures (2}, (3), (4). (5), and (6) demonstrate the results: obtained 





by integrating Eq. (35) numerically and comparing 2,, to the analytical 
approximations of Eqs. €36), (37), and. (38).. Figures (2).,.(3),.and=() 
were drawn for a depth of 20 meters and Figs.. (5) and (6) were drawn 


for a depth of 197 meters. Figures (2) and (5) demonstrate that 2, 


as a good approximation for 2, when & < 100 p.. However,. this: approxi-- 


Ie 


merion breaks down wien £ < @0 us. “This is because the effects of- 
fuerace tension, which are important in bubhles-of % < 30 u,.are 


neglected. Figures (3), (4), and (6) demonstrate that 2,,.. is a. good 


iz 


meproximatven for 2, whem £ > 400 1. Erom Fig.. (6),. 1% is. observed 


a 


that 2 is not a good approximation for 2&.. This- is: because. 2 


T3 cs 


neglects the effects of gas diffusion. 

The survival threshold for bubbles is different depending: on 
whether the bubbles were released at a depth 20 meters or 197 meters. 
faethe 20 meter case, bubblesvor (> 100 jj will®@reach the surtace.. 
while in the 197 meter case, bubbles of & > 1000 meters: will. surface. . 
This difference is caused by gas diffusion. heing dominate. over: the. 


decrease in hydrostatic pressure. Since gas diffusion is proportional 


to (P-P)/P, it will be constant for z < 30 meters and will only 


20 





decrease, for 2 > 30 meters.» faucsuubiies released at 197 meters are 
exposed to a high gas diffusion rate for practically all their travel,. 
while bubbles released at 20 meters are exposed to a. decreasing gas 
diffusion rate. 

Finally the bubble distribution function was calculated. for. each 
of these approximate solutions to Eqs... (22a) and (22b), (34) and‘(35),. 
Wo5), (37), amd (38). It wes convenient te introduce the imitzaie 
distribution function as the boundary condition. (2;z ,2_) instead_ 
of as a distributed source. Since s(z",2") is zero,. the solution.can 


be written as 


V 
ize y=” Bz 2 ) 1/(ae/an'). (39a) 
© 6 Vir O° Oo 


hor y = Vins Eqs. (34) and (35) must be integrated numerically to: yield 
f . 
‘ = { < = 

o, (232. k 0) Bore Vp ( < 100u) and f = 0 








; eo - 2 ‘D-z e 
— ; 
for v rot > 400u) 
| 1/2 
CO =Oe ae Ce a oz.) (39c) 
DPT GAG D-z | oo 


and if the effects of gas diffusion are not considered 


Use D-z ee 
"e poe ee : Ver eee es 


(39d) 








O 


heme V_. 1S computed using 


ws . v3 
tics 
e Dez 





O 


Preigee(26). 
The above distributions were calculated using an initial distribution 


Funct Len , (Zz sk _)> obtained from Ref. 4 .. The authors of Ref..4 
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determined the distribution of bubbles rising from "the highly organic,. 
anaerobic sediment on the bottom" of Saanich Inlet,. British Columbia. 
The depth of the water was 197 meters. The density of these rising 
ee was determined by using stationary echo-sounders operating 
at various frequencies in the kilohertz range... Figure 7a 1s a histogram 
that presents the presumed bubble distribution,,. o(z 2) at a.depth of 
197 meters as presented in Ref. 4. Figures (7b) thraugh (7g) demon— 
strate how this initial distribution changes as the bubbles rise: in 
height. It should be noted that height is measured from the ocean 
bottom upward to the surface. 

The different histograms at each height above the ocean floor shows 
how the bubble distribution is changing depending on the: model chosen 


to be used for the calculation of 2 In the model where gas: diffusion 


rr 
as well as hydrostatic pressure changes are considered,. the bubble 
distribution 0(23Z sk) shifts toward smaller size bubbles as. height 
increases. It can also be observed that the width of the intervals 


1s increasing with height. Finally all the bubbles disappear entirely 


at a height of about 100 meters. 


i 


In the small bubble approximation 4 @,(z3z ,2_),. where the drag 
a oo 


Very? the bubble distribution 


force is considered to be viscous and v = 
appears to spread out by shifting toward bubbles of bath smaller and 
larger radii. Bubbles with initial radius less than 900 u continually 
decrease in radius, while bubbles with an initial radius greater than 
900 pw will initially decrease in size but as the effects of gas 
diffusion reduce, these bubbles expand. In the turbulent drag model,. 


v= ¢ = o,(23z 52) shifts toward smaller radius bubbles because 


nie 
of the effect of gas diffusion and disappears at a height of about 125 


meters. 


meee 





When gas diffusion is neglected, the bubble distribution 
> = o,€23z 5,2) shifts toward bubbles of Larger radii.. This was to 
be expected because as hydrostatic pressure decreased,. the bubhles 
would increase in size. The bubble distribution calculated for 
(232 52) compares quite well with the experimental data of 
McCartney and Bary in Ref. 4. It appears as though gas diffusion 
was somehow not important in the bubbles they measured... These 
bubbles originated at the bottom, so perhaps they had some: kind of. 


coating that impeded gas diffusion across their surfaces.. 
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VI. CONCLUSION 


Calculating bubble distributions in the ocean by solving the 
bubble transport equation appeared to be quite feasible.. Although 
the solution of the transport equation can be quite involved,. a careful 
choice of approximations often permits approximate analytical solutions. 
Otherwise, the equations can usually be numerically Integrated... For 
the model chosen, analytical expressions were obtained for several 
approximations. All of these expressions yielded results that: were 
in general agreement with experimental data or with intuitive reasoning.. 
It was seen that the drag force model chosen and the inclusion of gas 
diffusion greatly influenced the results obtained... 

The assumption that the velocity of the bubble was equal. to the 
terminal velocity should be quite valid as the data presented on Fig.. 
1 demonstrated. With this approximation, ome can easily calculate the 
bubble Padwus dt anyedepth. THe Padius 15 Seremety sepemden > Online 
model chosen for the drag force, gas diffusion and surface tension... 
The choice of the drag model determines the equation used to calculate 
the terminal velocity. 

Finally the distribution function for a bottom source was calculated. 
The data generated by this theoretical approach agrees well with the 
limited experimental data available. The distribution calculated for 
the no gas diffusion model yields data that is in good agreement with 
the experimental data in Ref. 4. 

From the limited amount of experimental data avatlable,. it appears: 


that gas diffusion is unimportant in bubbles originating in bottom 
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sediment and transported to the surface. It is felt that the type of 
gas in the bubbles would greatly influence the amount of diffusion that 
would take place [Ref. 12]. Since a fluid, super-saturated with the 
same gas that is within the bubble, is the only fluid in which. gas: 
diffusion is not considered the dominating effect,. there must he some 
mechanism inhibiting gas diffusion from sediment bubbles... Further 


investigation into this phenomena seems warranted.. 


25 








[ Sahist 3 


SNIGVY S31IEGsNd TWILIN! /WOLLOS 330 LHOISH 
OV OF O¢g Ol 


D3FS/SYUSLIW 0'0 -=°A 

SYSLAW Og -=°2 
ALIOONISA IWNINYSL OL HOVOUddY SHL 
‘avy 31gEnd sO SLINN NI 
GIIBAVYL SONVLSIG SO NOILONNS V 
Sv SS18gNd SHL SO ALIOONISA SHL 


uw) 
N 


"@ 
T@) 
ALIOO TSA TWNINYSL/SALIOOTSA 





SO <a oe ee eee ——= — —+100 


25 





gZ vinstj 


SYHILIW NI 
WOLLOE 4d0 LHOISH Od! | 80 


O 
O 
aN) 


= 
¢ 
(™ SNOUDIW NI Sniava as1gan¢ 


NOILA NOS “WOlY SWAN 


LOVX3 OL NOILMA IOS 
WOILAIVNY SIVWIXOUddV 
JO NOSIYVNOO V 


SYSL3SWOd-="2Z WOYS N3SIY 
SVH 3188Nd SONVLSICG 

JO NOILONNS V SV 

SMIGVY 3188Nd NI SONVHO 


O 
O 
O 


008 





“415 
ee | 
ee O00! 








i oe 


SYSLAN NI WOLLOG 350 LHOISH 
OO? OS! O00! | 


NOILATOS WOISWAN LOVXS OL NOILATOS 
WOILAIVNY SLVWIXOUddVY SO NOSINVdIWOO V 
SYSALSWOd= °2 WOYS N3ASIY SVH S1EdgN SONVLSIG 
JO NOILONAS V SV SNIGVY 3188Ng NI SONVHO 








10G2 
Y0SZ__- -_—__- - oe 
NE a eee eS 
= 
~~ 00k ec. as 
———— eee 
a TOS ae 
_-—HO0R 





O00! 


NOYOIN NI SNIGVYE 318EsNnge 


OO0S 


(77) S 


OO0O0V 


OOOS 


28 





4 OANIT I 


SYSLAW NI WOLLOd 430 LHOISH | 


002 O¢| OOl O'S 
a 
ee 
7008 2 = 
= 0s a 
-—— MOG 








“7 708 tee : 217 
ee ae -— —————-— =] 

77000! 2s | => 

Y000! = oe | ‘NOILNIOS WOMWSWAN LOVX3S OL NOLLMT1OS 
an TWOLLATWNY BLVWIXOUddY JO NOSIWdOD Vv 

Z SYSLSWOZ-=°Z WOYS NSSIY SVH S1agng JONVLSIA 


“a JO NOILONAS V SV SNIGVY 31esNg NI SJONVHO 


JOOS 


OO0L 


NI SNIQVY 3 1ssNe 


O 006 


0002 


(77 SNOMDIIN 


OOE) 





300 


BUBBLE RADIUS IN MICRONS (2) 


50 





lO 
LO 20 50 4.0 9.0 
HEIGHT OFF BOTTOM IN METERS 


CHANGE IN BUBBLE RADIUS AS A FUNCTION OF DISTANCE 
BUBBLE HAS RISEN FROM Z=-!97 METERS A COMPARISON 
OF APPROXIMATE ANALYTICAL SOLUTION TO EXACT 
NUMERICAL SOLUTION. 


Faure. 





5000 l= > ea ea 


BUBBLE RADIUS IN MICRONS (4) 





lOO 
20 SO lOO ISO 200 
HEIGHT OFF BOTTOM IN METERS 
CHANGE IN BUBBLE RADIUS AS A FUNCTION OF DISTANCE 
BUBBLE HAS RISEN FROM Z.=-197 METERS A COMPARISON 
OF APPROXIMATE ANALYTICAL SOLUTION TO EXACT 
NUMERICAL SOLUTION. 


Fisure 6 


G 
t- 





DISTRIBUTION CF BUBBLES ORIGINATING IN 
SEDIMENT AT DEPTH I97 METERS. 


NUMBER OF BUBBLES 


400 600 S0Q 
BUBBLE RADIUS (.) 


INITIAL BUBBLE DISTRIBUTION $&(Z,2) (REF. 4) 


Figure 7a 








SA 1asnd 40 YSEWNN 


BUBBLE RADIUS (/4) 


Sa 1Esnd 4O YSEWNN 





BUBBLE RADIUS (4) 


iN 
tol 





=_ 
= me ee 
a» o oe cate 


3 
BUBBLE DISTRIBUTION AT HEIGHT 


CO METERS. 


Figure 7b 








S3188N¢e 3O YsASWOAN 


500 700 
BUBBLE RADIUS (44) 


300 


lOO 





300 900 (00 


lOO 


BUBBLE RADIUS (AL) 


BUBBLE DISTRIBUTION AT HEIGHT 


SOME TERS. 


Figure 7c 


3H 


NUMBER OF BUBBLES 





d i po = , 
[50 S30 oO) (0 Se: 
BUBBLE RADIUS (AL) 


NUMBER OF BUBBLES 





BUBBLE RADIUS (4) 
BUBBLE DISTRIBUTION AT HEIGHT = 75 METERS. 


eee SS Se 


PI fel fel Ie 


WN 


Figure 7d 


35 








a Oo 
Sa7eEsne 40 YSEWNN 


600 800. 10.00. 
BUBBLE RADIUS (2) 


400 


200 


Cj om 
sd 18snd JO YASWNN 





BUBBLE RADIUS (AL) 


IOOMETERS. 


BUBBLE DISTRIBUTION AT HEIGHT 


ure /e 


~ 


Fig 


(O 


cr) 





NUMBER OF BUBBLES 


NUMBER OF BUBBLES 





200 400 COO 800 IOQQO 
BUBBLE RADIUS (AL) 


BUBBLE DISTRIBUTION AT HEIGHT = I50 METERS. 


Fisures7 © 


400 800 {ZOO ISQO 2000 
BUBBLE RADIUS (x) 


BUBBLE DISTRIBUTION AT HEIGHT = I97 METERS. 


Figure 77 


ay 


——aee gee en 








COMPUTER QUTPUT 


oe O00 METERS. 


OFF THE BOTTOM 


HEDGHT 


— DMRODOAMMMOMoOhKAOMMMM++ 
ry NANMO SAM SIMMS AN MMMM Mee 


Fee vss cf ee 0 a ke 
BOOSS DOONAN MM MOCOOGOS 
MOO ea 


LAL HORM OFNOR MME MAaOoNOtONn 
MOON WORALAYEMAMODIMMhDOSFEKEKO 
FLA MON EDMANCOMINOMAAMROSDOWM 
~_—) ®#seoeotiegeeebeouoewoervsieev et egeoeeeoese és 
LM IANO HAA MOUDO HON MOWINA COLIN 

ADDO ANMOANOGMSONNODAMOMOM + 

SPP ALAS LAU EEN 0D WOOD DOE EES OOO 


FIONDOONNUNANAS ON CUTS MOM 
INIA AI ODD OIA DOA CIA Mmmm mae 
MOO ft EUENLALINA- ODE EE ADDO DHHO 
~Teeses#senreteetrteeteeese eee © @ @ 8 
CLOVDODQVONAANMNNMAQAMOODOD00 

MMO OM eH etes 


WY NONMO ODOAAIO SAD ONNM™ASPNASwO 
D VDOOMAORANROMMNAIST OIVNODDOMANSS 
m= OODAHINO DUA ONMDOFAMOODMYO tev 
QOeeeeetrteeereereeeteeree*eee# @ 
GOVDDOFAOTHKAMOMTRHALTEHRE QAON 
XK FTN DOOANFODANFONMEOHODO 

AI AIAIN 69 0909 OT tS SF LILI LN WO WOO OD 


QODVQVOOOD OO VOeCeOoe Gee 224 
QoO°COCO COO OSee ea eOOOoe eae. 
ODVNDDAVDVDIOOVO AVIV OVI OVIOOS3NVOO 
wi @ ®@@@8@@e@e © 8@8HhmUhOMmUcMHhCMFChCUhHMUh/} HO HF FO 8H HG Oe @ 
KONN NOD BUNDNNNNDNOOMOONO® 
DOMONDNODAAMTNOEADOMOM-MOMSEYOR 
PPP PIN MIDI LN OOO OM EE CO 


C000 CO FHA HOODOO HAH HOOD MOMWNWH 
NN AMMO OO MODS + FLNINMA MM MAM MO 


Beats ck at ae ME) Moe e 
Cec oe ER ASSeOoS 
MOTO Aas 


OTDAINDOADNVOMNAAQ\NOMmMthyod 
NORMOEKOMOUOMTLAD SFANMOD+tOHO 
HEN APTMOOOMANADMOANMMOAKOON+TOS 
~J) @e eee eoeeoveeegeeeareeeegee8 & @ 6 
LOr-OMOMNNAP NVA DDOHONNAN 

BIHAN MSTNOODNHAMODOMMNANMBROS 

CIN OF ET tt SUN LALO UN OD OP SS 


Gr OY OS UN LOVED SP PERE EE EDD OOD D 
ro PE Pe LOLI OD OO PN et EE IEE 
bt C9 CO TD DM Hh SACO LO OAD DUN IN OO WOOO O 
T F®eerteereeseeee eeereeet®t_etete Ge @ 
AOOVODDOOAAGFEEMANNODOOOOO 

NNN 


MOPOPFOODKHROMNAMEM ODNDPOHONMOO 
HON HONG DOMSFORSeAHMADHONY 
mOUSNDNONLFUISEHNUOONDMOAAMOAM 
ws © © 8G eer*eeevreereef#eeeereeer#»ev#ee @ 
LOATH MODNAMAMDODOOLOSANAOD 

MOAFINDEHODAMSFOURDDAOGOSAPIDADAG 

OOS SP SS FUT LAL NLD UNO OO LORIE EE OOD 


eel ele elalel@ la la Sle) Sle Se) @) ai ela) aie is) 
OOOO GOO DO OOOO DOO DTOO0CO 
QOOOOODOOC OOOO OOOO ODOO0000 
wi ®e@ @@@¢#8#@¢@e © eee eeee?e © ¢ © © ¢ @ 
KONWMN MD ODNOWMNNN OM VDONHNODOMOS 
DOMODO SAMPUIMS DOMDOM-OMSMORF 
wt PP SP PLANNIN AEA OO OE YM COCO 


38 


ine 





SOROCO MERER'S 


HEPGHT OFF THES80rTon 


LA OSFM DAFT MOAMAMAANMMMS 
TM AWNIDADMAM DOAN ST FOOD 0000 
—ODDORNADRNAAHKNOOS $+ DDDDDDO 
Teeevsvleesveeere#erete tee? 8 8 © © @ @ 
BODDONNNDAGDFSEEMMOMOHOODOOOO 
NNN Aa 


MON AMDSMHAMDHOMNHNODMOMO O 
MADANODNHNADAEOSFTNOTONY YAMA 
Mm ODDO MNN AM SHOWIIATAMM OOF FON 
—w © @e ®e@ © 8 @e@ © @e © ee fF 8 @ ee © © @ & @ 
LOROAMODTA-DONEM OMTMORDM cw 

MON MO SFUAD DW HAOMUIAD AUC Ort hc 

FLA LA LALA LIN LALA LS UNO OO WO WO WO PF FF 0 CO CO cD 


OMOMODAWNMODOr 
AODAMANAIIAIN OO 
—OODDDDVCDDIVIIONTYNSSTODDDDOND 
Teervreoesvoereoeetteoe® eer eeeeet @ 


ALODDCDVDDOVOVVO OF TNNMODOOOC0OO 
aod 


—Y LAWN OO SF ON OM UO OO 
~ Ott DONVADANND 
—OODVOVOD DOD OMALFNAONM OWA 
ODOeteteoeverteretoeevs®*eoereee et © eo © © © & @ 
GOO VDOVVDVQVIOIO AO ODNDTFODDAHNONGD 
of DDODSOAAAMALOM 

AAI OO SF TULA 


ODVDDOOOC CO OCOCeOoascOoOeaecoc 
DQOODDVDDVOOVODVIOOVOVIODVIOOVTOO 
OdQD0OOOCOOOOCGCOOOOCe@eaeooe 
—_ © © @@e@eetrt feet ee ee © © © © © 8 8 @ 
KONWMWNO DONDNNMNNONDONDODNOSD 
DOMAOAQKDOSAMTOAM—MDOMOMOMSMOR 
SSF LALA LAA LALA LA WO SO WO Be Mf 0 0 


AADODOOOOMNUNAT ttf aad 
N PE DADDAARAHOOODOOREMY EP 
—SOODDOVDOOVOODVOOYT SS DDDDDODO 
Tecsocv7wese3v2+#v#eereeteeerteeeeee © © © © © @ 
BAODVDDVOOIDNMNAAHETECMAMMODVOVDOVOO 
NNN AAA 


OtAODANEMSNO-Or-ODMSTOAEN 
TMYODPAOLTSFAANRDOMNODNDAOR 
SNE ae cee on ae ees 
—_+ ®©e8 fe © eet eeteeeet eee ee @ 
MODDAMNWDWOANDAMHMOSTOA-PANK-EADN DO 
PF KHATFOADANODONHAAIDIDANDO 
AAAI AI AICO MOO SF AE PULA SO 


COTM QO Qvh MM Ferg OODOMO® 
Aa DDNNODODDDVAAKAMNYIMNANMOMUIW 
mO-P MM OM OMMEMODOOOO-RY--mO 


Leeeeere#e+reteere eet eeereeeeee 
AO ODONNWWWALIANANOIMOOMODQODOO0O 
NNN A ac 


MTOM AIST DRACO NAD OM ACUO MIND O co 
SSO ANA PAM QI DSL ONINO NS 
LIMON OMNIM AMM SOM ODN MO FUANS 
ws @©@e eee eeere ®t eet Geereeetet @ @ @ @ 
LNDDAMNLTODTREAMANDNAMODONS SO 

NDON MAOH AAMOODATHRAONOMATY 

A109 FS ES LALA LAN WO SO BO PY OO CO 


OOO @C@O0GOQO0000070O00CO0Q0000 
OVDDVDV DIO ODOVIO OO OVO OCO0000 
OOCOOO0O VEG OC2OCoOOO000090900 
| 
KOWMMNQDOOMODNWNANONMODONOONOO 
OOM DDOAAMTAMSHAOMOEOMSEH OR 
SPATE PT LALTIALALALA LAIN OO OE I SM 00 


39 





15.000 METERS 


HECHT GFF THE BOiimEn 


WOMOPADNON DOAN AOD OM™ OO 
MO eEe~MSESTTKONMS DOOMMMADDDAOMD 
mEOM MRE PE STINO NNO AIEEE EEE O 
Teevrec5ueo5usoese3#es3eee0eeenstee 8 © © @ @ 
AODDONNND ODUM MAMNAIADODOOSDOO 
NNNAaA 


ADOFMADDNNOONASREONSHE NO 
MA ANMOIM OOANINDO MM OTVOMHAMASTO 
me MAIMO FNM KY HoH OONOK;&FORhkeOtTO 
—_J ee ®@ ee Fees eeeereee#ee © 86 © @ @ 
SCTFOAMKHVONCDMNMYFOAMmMM ODO DOMNODA 

OFUVO DO OMDAMTFODDYHOMAPNOAS 

PULAU UV UNO 50 0 00 OP ME ™ OD BOO'D O 


meODDOODOCVOOVVODOVOCVOCOOFLE SFr Oo 


TL eeevrteeeee test eeeteeeeetee ¢e 


o OOS OO0OOOCOOOCeOOOOCO0O0CO0 0 


WY MAOMS + 
=) LUVAINS OO 
—OOVOOCVOOCVOOCOSG CGWe@toa en: 
Oereteeereeeteeenreet © © & ew © © @ © & @ 
GLOQVDVDOVO OVO OVIADIOVOVOABODIFMWAWIDY 
oe PT NRA 

MOUND + 


COCO OOCOCOO COGS Oa oe eeo 
CQOoVVOOO OCU OOCSe@ee Oo eae2za® 
COOCCOCOOO COOCOOCOCOCee eee e sa 
—_— © @e@6e¢'®e @ @eeee8e @ ee eee ¢ © @ @ @ 
KOMWNO DAWMNOMNNNOWO ODNOANA® 
ODOM ODO HHMSPIMHOMOM OMS Or 


UVF PTHHDHHAD 


N I= D Os OY O Ov O08 OVO 
QO ODDODVDDOOCOC ON AAA ye EYEE EEO 
Teeevoevvreeeeer eevee eeveeeer eevee 


AOVOOV OVO GVW ONMNMAMNMYNDOODN000 
— 


LIANILA DDO DDMOAA! 

Soper sts AN Ost 

mo ec oemmmoOo OKI nee oon. 
—_ ©? @®@e2 @@ 8 @ ee © @ ee © fF @ © © @ 6 @ 
COOWDOD ODD VDIAVDOONANMNANS FOODOO 
APIAMOWN OAT ce 

ANAL OOO SS 


Om OS tad HAAtHOWDDMOwOo 
A MOAMOO OOO OF FPOMNOMNYAMMAMNAY 
m= OOOO TT STOODODDONUNDODODODOODOO 
Teeevrsevseeereee Feet erteeeree @e ¢@ 
AODODOST TS TAARHOOCONNNOOCOOCOOCOOCOCO 
NAN AAaS 


Were CAD SOUS AOD ANS ™ WMO OOM 
a OOTFTONMODODANDONNESADODMNN 
FOROOCTOMMAMSFNTNMODO0OLTOON 
— @® @e Ge e@eeeee Ge eeeers © @ @ © © © @ ¢@ 
CLORODOANAND™~ AMATO NEAMon 
MDVD ANTFNDONODOMDONDAIASTAD 
CVC 0 OY SFT PF SUL UN 00 OD IY OOOO 


OCS Bae CeQOQGQOQaQaoQaaqagdaQaQIooO 
OOOO G@e oO ee eee eeoOOOOU0O 
OODWeO OOO Ge COU OOOUeae@OODCOM 
—_ 2@ee%s286860808#eve#ee#h#f8f 888880666 @ © @ @ 
KONWMNOVDONOWNMNWNIWNODONWDWONOOD 
DOM DHDOWAAMTNEHDOMOMONnTHOre 
est Sh PUA UU UA UV 0 0 OO OP i COO 


40 





100.000 METERS 


HEIGHT GEF THEMBOTian 


QANDAHOHAPOMNODMME NOS PUNO 
M DAHDHPWNAHMOFVDRODOOODO0O 
m4 ODDO D ODDLY CODDOPr DOREY MY [~O 
Teeee7s3+ooreer+3+oese#eeteteeeesee? et eee @ 
AOOOOSFHEIHMOMOAMAHAANNNOOOO0O00 
NNN A eted 


AI SPUYV UVLO 0 OOM DDHROOANMMSNW 
MADONNA AHMODD+ODAHMOSPNOOrND 
et tOMMOST MOWPMONMOT CODON 
—7 © @eeeeeeeeeee#eeeeeee#eeee? ¢6 
CLORATOMNNAHODRHOAHEANSORNNOD 

RMDOANM POM BASF POANMMEOND GO 

PUIVUVOODDODOOOPEE EOD DAAARGHO 

rod 


© 
0 @) 
MOND ODNANVODVDAVDDVODVONVO0O00 OO 


Tl eeoeeerteeerereeerteee3++reereee tt @e tee 


LAODOVOVODOGOV AOI OVOOO90000 


Y COL 

a ~O 

—“ODODODDVOVIOOV DOO 0OCOO CO OW) 

Oeeerkereeetegeeeeeeeee7etf ee ete e 

AOVVOVOVOIOOVIOOO DOO ONOO OOOO VOT 

ce ae 
oH 


ODVOO OOOO OGO460460040000 Gee@ 
OOOCOOODO OOO O2OOOOe OOO Cee 
OCOOoOMQWOOVWOODOOO Cee Coo cee e2e 
—_ * @eeeeeee#e#t eee @ © @¢ &¢ © @ & & @ @ 
KOWWMNNODONONWNWNMNONMODONDONHDO 
ODOM DHOHAMEFUOM RO MOEOMSHh OH 
Pt tet SUID UIA OOO OD PY i OO 


DADRA 
N me Hed ed 
aOODWDODIOODVOVDIOVVDOVOKEEESPE EO 
Teeeevroeevueeneveeeeseeeenret tee e 


eelelel@eleleleleleleleolelelelelelelelelelelele) 


DOAYQODLAD 

N AaINHOON 
OOD DOODOVDOOODV VAIO OOMN Os feip- 
wi @@®@eeeteeereeeteteertereeee ee & © @e © @ 
CODDDOVVDDOVDADODOOOO LE EeTNH 
LAYS ONS 

a ONO (} 


P= PhP CO CO COWL Mh OO ORY Er 
Ste eH I OOD O00 OHH eH. 
mt OWL 0 0 DUAL IS BM OOO LNMLALTUN LN © 
Teeeevrereeveetrtereeteeeterereeeee 
LAODDONMNMEYEYDOWDAINNDOOOO0OO 
md md ed 


KHODDAMSFOVINDMNAOMADENOONA® 
a OM PRON QAOMU eH AHOWNNOMAITY 
ROOM MORN OANWITWN Dam OS ES SOOO 
—4*?eeeeeeeeet*eeeteeereeeeeee i.e 
CLODMWMNAHOMAHOMMF DOOM OANODMTOW 
ADMIN DANMORN DONMAMANSAMS 
QIN ADV OSS PUVA OO OEE HO DODD 


OOOQOOOG00000Q00000090000 
OQOQQOOoOeoQgoOeoee QQOOQDO0O0000 
O@G@20O0 DOO0OO0 CVO OOO 00 7000 
—w4 e@eeeee65oueeecs85§5eeeeeesbknRneerteeeee e 
KONMNNODVDONIDNNNWNDNHNOONOOWNDS 
DOM DDO HAHNENEHROOORPOMT™M Ort 
SPP PLA ALU 0 OD OO I EH COO 


41 








1253090 METERS 


Orr Tht COC ToD 


HEIGHT 


OM Wear com FFODOOMNDORNMST Ss 
M DAADHHONVOVUDONN VARA AAAS 
mt QIN A AOS PFI ON OS SOO OOO 


Teevseeee7eseeteeeteee ee @ © © @ 
ADODOVOSSTHFOOCVWVODOONNNOOO0O000 
NUNS acted 


MO FOO PFODNOMMNMMNAHOSFOUOMS 
TOS FO DDOOAIAIIM DANO NOM OM OF OO 
HORM MAHOU OOODADOONDONNAD 
—_/ @ee@®@@eeeet @eee0e © ee Ge © © @ @ 
SMFOHOO-O FHAANANODOFONDMOEUIAA 

FOF OODRHON FND AHO DANHONDOO 

WOODDOOEPMEE Fh ODDDMADNAOOOCeR 

asd od 


LO LOLOL OOo el ela Slelolelelelelolelelelele 


Teeeeeeeteteoeereeeeee tee ee @ 


COO DVO OVMODVODODO9O9O9000000000O 


VY) 
= 
MODODVOVODVVOIOD OOO VDVOVOOVOOCNCD 
Oeeee#eeeeeeeeeeeeeereeee¢#e¢ @ 


LOCO ODOCOOD0O OO 22022 oOe2eeeo 
ao 


OOOO DOOOCOCO OOOO O0 Gee eeaage 
DOOOQOOCOCEOC6G400 CeCe eeooe® 
QOVDOOCDOOOCOCSeOCeeoOooeaeeoae 
—_ ®*eteteeeeeeeerkee eeoeeeeee e 
KOWMMWMAVOONONAWWAO NOONOONODO 
OOMDDNORAMSFOMHDOMO-EPOMShOR 
Pe SP SPA ALAA VAL) 0 OO WO OE I COO 


oO 
NJ N 
—aODdDOODVDVDVVOVWIOVOOOVOOOVO 
Tececeoeteee5ueoeeepmeweeeeteeeteeese @ @ @ 


0 OOO OIG Cae 20 COGS @QOOCOOO00O 


COON 
con 
Molelolelaleletoteloletoleloterelelerotetoron 
eeeeeeveeeoereeteeet ee ge & @ 


zal! 
CNS SNe ae Uae Oa 


— rm 


WD EA CACO OD DO FFF UEVUVUN OO 0 00 0 
| DRM ODO MOO DODODAHAAAG 
mt OD O00) Be FPR PE LALA LIA CO MOI OM CO CO CQ © 
Teeetereerereee*®*eeree tee eee ee 
LODDONNANNMMAMODODRAAOODOOO0O 
aoe 


NU FODMAMOAOFODOMAMMMH DAO 

AOD FIV AIFS OSF ODM (ODOAIS 

KH COOND DONO TARDOATFOODDO 
J eeeertevoeveereeveeeeveeeerevoet eet @ 8 
CODOMAS NORA DNODON MD HEIN es 
MOS DAHON OMNFHODOSFHODON 

NCI OO PLEX UAL OOM i OO DDO © 

4 


OOO V OOOO OOD OOOODO OOOO O000O 
OOOO GO0 EOE O2OoeeO00O7TO00000 
SGO@Geaeo0IGogeooeeQoooog@goo00o00O 
—71 @©@e@eeee8e ef @eee @¢@e @¢ @ © @ © © @ @ 
KODNWMWO DONO NMNNNOUMOOMOOMOO 
DOM DDO AAO PFUIMDOMOMONSMOFS4 
SSS SP LA AVE UAV DD) UV UN) 80 50 0 BE I OO 


42 


150,000 METERS 


OFF. THESSOTTeM 


HEIGHT 


DANN ADHAAORMANNAS SPOOR DD 

MY HADTFNNOMNTMNDOOOHAOOOD0OSCOD 

MOST PPODAOFTNOOOMMNNUWIING 

Treempetreerseeeeeeweereeertee tee ee @ 

LOODDOMMNMNDODADDANNNOODOOODOO 
aos 


LAS OM MOF FOVUMAMUIENOADODTNMOOeR 
CVO FAM FHOMAHONOUIVIAKE Sheol 
COM APHID ODMODMDNND S$ DOOFN.00 
+ ®*eeeeeeesetee te Ff 8 © Fe © Co 
Cr TFFOWMMD YALE MOOT NOFHMAIRAAW 

SrICTV FIND OD AMN DO ANB ANDAM OOS 

OPE ERE EMH OCMODAONNDNHPNOOOAAANN 

Saeed 


—QOOOOC0 SO COO OOOO COea 2 OGCe 


TeosveeevrsetFewvreersee ®t eree e 


OAOOO00S8O OOOO OOCOOS8oOeoQo 


WY 
a) 
—OOOG0OOOOCOCOGOD2OO0O02O002e@ 
QDQeertete FP eevee eer ee ®t 8 eet eo ®@ 
TIOOVDOOOCVO@VOOOGCeeQOOeCeeae@ 
~ 


COO COCO OCDOCV EO Oe ee eae] 
OOdO OO OOOO 0OC COO eaeoe2eageo 
COO OODOOOO {2OOCeOOeaCewo 7a s 
—wJe® ® ets ® ee®eereee epee eerste @ 
KOWMNNOODOBDNDNNNNONODNOOMOG® 
ODOM DDO AAA AMAOMOMOMSTHOR 
we SPE LALA LAV UV LAL UNO WO DO BH BO 00 


N 
SOODdDDODODOIODOVODONVO DOO 990200000 
TL eseceeeeeeeteertee test ee & @ 


LDOODODOOCOVCOVODOCOOO 0000000 


N 
J ©e®eeeeseseet eet © Ge © eee & @ 


COGO0GCG 200 OOOC COO OO O06 OO000O 


COD AINN SARA ARANDA 
wont NN ODO DDNDOOOME EY YE 
RFOMODOAIAA AMMO OO SHA ANNNNNAYO 
Aes Spugtiaeiata nat 8 item en ces ”  f eka? 


DODO DODAADANSPSESAHAAOOV9DOOO 


PIAA OO a SAODMIARAMONO ENS 

ALA FAD SOAR MADN ONO AO 

- DOS OOMON SERNA DONO LON 
—_w7 ©ee © ee? tvee et ee fF Ge Fe &€ FG @ @ 
CODD DODNFOONDSEMIVMA AVE MOOS 
COMODOSPMOHDOASTHSTOAONS 

AICO SF FINO OP MOO DWOONAe 

a | ed 


OBeaQaoeeo@oQoQooQQoodvdoQoOdO0O000 
8 (Se (SS (ae SoS a el@elele lel @elelele @le 
OOOOOOO0CO2E2Oo9 Gea oOOO0000 
—7 © eee eet essere eoee#et ef @ ee ¢ © © © 6 @ @ 
KOWMNWMNODOWMNONNWNOINDOMNDOMNEO 
DOM ODDAAMEINMEDOMOMOMSROR 
at SS LALLA LDV LAL LIV UN OO WO PE CO 


U3 





R75..000 METERS 


— 
7 


HEIGHT OFF THE BOTTOM 


OHANONMOK Ye t$OOORAMINOON Sd 
MH OO0MSHONDTOODDOOAREEAERE, 
ORM MUNN US t Node fF Pe) 969095. 
15 Se eeeeete @ 8 ee @ @ @ @ 
SSOOONN ANNO OSIAHASad0065 


tODAMOMOOKAFOVUONNTODN DOR 

MO APODMOMAIDBN OM FINE NOOO AQO+S 
EMOFODO OW AWA AMODNAEMMONAHDOON 
a *+¢6t+te60e866 ea esevetbevneeeesevnre6s 
CINODNMONAODEMAFHOKMDOMODYS OO 
FODDIAMMAD AMM OFONSHOKSTHAO 

DWONDANDHADDWVODO SHAAN N GMI se Sewn 

Fm hy ed Sd meh a et et 


RINOO CO OCSOOCGE2 GC OOCOoO ©0000 


TFeetovoteoeveet 6 @ ee 6 @ 8 © 8 6 6 


AVD@AOCO @QQOOCVOOSSOO0WV@0VGV009R8 


VY) 
<a 
HMHOOODVIOOODVDOVOOVOOOOCDOCCOO0OO 
Ci eee ee eteveeeveereeteeees © » 6 
pe Neral SS 


OOO 37009200 CeoOCOO0os eGolTteae 
DOOD IOOVD0OOOE00C OGG eo2oOo 
OVOOOVOSOSOQ9DBDOOS99O9N0OG000900 
er ee ee er rr 2 2 ee | 
KOINNWNODOWNMODAWUNOMDOINAONAGD 
BOM DODO AAA SAM AOMOMOM Sh Oe 
SP eS SLO LA LSU LU OOO OE EEE OC 


ieee 
MODDODODODDOODOOCOOCOSOO0C00090 
ce ee eee ee eee ee ewe ee oe 8 


OSGOOD DSDOCOODOODCO0COSCOOS 


N 

KF-ODDDODODOODOODOCD0O000000 
a ¢é et vd,seebCeboev steve see bid 6 6 8 
Pease OoeOOCoUOOoOOCOOCOCCOoOSC 


QIAN DODD DWDWMDMONINWNAIAW 
rH COD LALAIN NIN Mama AINA AULA 
= OIODDOIONNND 00 000 AMAA HRHO 
Teer esvsssveereeeeeerteseeeteteese@e ¢ 


A OOQOIMANUNNNNOODOOD00O00 


DODtTOIDVONOOKOSOHNOS 
OMmAODODODEFNODODAMNWNMwaAaSt 
eS OOan CA aim Com mom 
wJ @e@@eeteeeeeeesvseeee tee ee 8 
CLOODSODOMA MALTON M LPHOdmaANnGOOlaet 
GO COPA SAID UIA SPA PMOANOW ED 
AKATIVODOMH OD HNOOKRAN OCD 

ped df ed md ed et md 


OCS ODD 2OGOOO OOD OODODOM7O 
QO .1O OSGeo DOO IO OOOCDG MoO OOS 
OD 100 70000200) O000000O008O 
wi @eeeeseeteteeeetee @e 8 © @ @ © 6 @ 
KOMNANIDONO NINA NOMNDOMOONOO 
OOM ODO Sel SLAM DO MOM OM FY O et 
St th FINALL 0 00 OE EEE OO 


uy 








LO TS20Ce> Meters: 


— 
— 


OFF THE BCTTOR 


HEIGHT 


ORE FO MAM OMOrMMmaAON Ss OMnM 

DANOOHLEONNMNAHNDODIOOO 

Am Hy} 4 Qe SOAS 
ee ee ee (©. @ac?. Sa"G © 0 8 SLO 


@ ¢ 8; 4 ena } j 
DOOD HAHADODMMMOSGOCOSGSOOG 


DADOPAID MOOD Oatask FR MNINDMLAO DO 
COD AM OO AT QE DOCOMO MOM Mah 
IIE ASO SOMO Om Mm FIN RdOM ODULNAO© 
-jebeoveveteeeov sete bei gEvte ti bw ES 
MOOTANNOSNOMcCO Mh MOIDiAt ahs Ot 
DEOL OANAUPE ASSAM MOMNAN 
COLO OTT OEE LAIN OUR OO MDOORNHN 
ited Sea et ted et Hedi HOI ROAINAIAIAY 


I PLe Salelelelela 2 alee ls) le (@lelelnle ol. 
Téesetbtepese_beuvdseoeeovoespeeeoesb_rueasaeeed 


pMeol@lelelolelelole(<lale)—~lelelelelels< lele lois) 


Y) 

eg 

NO OO 8 ae ae 
Meeevuxv4e54uxete eeerteeteee 


I30690090060050090000000 


eelelslolalelelel@lelelelel@lelelSlolelsiele 
AODVMOOVDVOOCOODWO VCO Vee eaea] 
Ol S 1 Sle 1 @ 1S el P lel el elele ea (m/e lele lm) o1@) ole 
w—wi e@eeeeeeveveeeteeeeeert* eeeee ee 
KONMNINODOD ANNO NWO NODINDOMNO®D 
DOM ODD Het gPareDOMOMHOMSHMOr- 
PPP PUA UVLA UAVS 9090 0 OO EM M00 


NJ e yp fe 
HOBO eOOoOOCOSoOOSGOOOOCOOO 
ee @eeeeee##e se ee «© @© © @ @ @ § 


F one doc see daeodcdaad 


Poadc00e99c0000000C90000 
jSeeeCeoeeeveevneveeeveosvredce & § 


Melelelololelololelelolololelolelalelolotelole 


COC ST EN HOO OO LOO 
_ aod at NAlOvon gist thd tt 
MHOMBWODOOMVODOI ODDO RAHODOOOOO 
se ee 2 ee 2 2 | 
Ole l®lelelelel Ll e1S1Slelelelelelelelelele 


MOM DAR AIOISIS DM CMONDQININ 

PD ODO MAMA (IVKM FAI DOW 

= IOS doows LA COLA ST LAP OF ONIN OO 
— J ®@e © 8 @@eeeeeet%e68 &®&e @ @ @ @ @ @ @ 
LODO OAD AP RAAAM AITO ODOMNMODN 
NADWDEODMAAMNDODAAAY 

DN OM DDAAIAIM SUA EE ODO 

Med HA a Ha HHA 


COOVDQMOOCOOCOCOCOOO9D9N00000°8 
ONO DODO DBDOOOOVOAA9090000 
DOOCODVDVOVOOOVOAMOM9DO2D92000O0 
—_/ © ee ® @e © ee ®@ e © © ee Fee © ee 6 oe 
KROWWMNODOMOUNWNAWI ONS DWI OWM OG 
ODM DHOOMAM PUES DOMOM WOT H Oe 
PPPS SUOMI NU. 0 ORES OC 


US 





COMPUTER PROGRAM 


~ 
— — 
QO nana 
eo wd eed 
~ - “ND 
Oo WW Is 
N NWN w+ 


OLN UY) eet wpe C1) ee 
am OO wuld | we st 
Satu tm 
OO ett ¢— 


Foe NG AN 
oe eM NTE 
SOY & & > 
~OIYIL DO 
tm ONE 
mm OJ me Ue 
Nar ODTEN 
} PHAR N 
Dw etn 
eer ULV OMIN 
THR NN~I> 
} TO ers IS 
<b me ee 


wee LT dt 


aN Cran 
DA~ Ca he 
Ht #ONM>nN + 
JO) o~ -CO——_ 
fe Ko tran OY 
tUxXNOININ> 
te ered rt pw Ht 


Nee 

= Lu 

st pat © I we co 

rs po ® at i ~~ one 

we «if th eZ “t © | m © ~~ Y 
~~ oC : Us oO mm WY i] e fe a8) _ > 
= _@ O om }K- ees N “ _ <I 
~~ WwW OQ ah ru me oO oy WY N So eC 
sl LL =) alas & owl Lid eas re AI ow 
~ =z hm . oe ue (Lhe at a oO Lu Re <I 
=— “se <t ¢ ud O = .O , = =< ei ot a > 

=m i we Ww < ™ ke bk cae ~ be ea +e, uu 
mam «a <«l Cac ea o< 4 Oem ny WO -. < a ~ © e) 
Se x J LL} - Q Zee Win, 4 a ~~ Oe A <1 
>= N Jue ote - | | Co 7 - if > = ce 
% SDS) = Y oc rt lu Ae Va iS Zoe NN ~ & > Ss 
N ~~ <A mY) © O ce Ge tL ee er > G&G Be ors. ~ b= 
NO ht U. ~ - Q Me & i - - O ~ Y 
~~ LY a Us ca ~ Lhe rH = — Lt} ey -» OQ —~ Oo 

| “kbieeel OF wa Lu it t aaa <tt- tl oe rw tu J ole MJ 
= 2 ml Cr D | e FY eG ages —~cy e@¢ I << x Ww 4 a 
<< D> Be = 7) i Oo 4a oh Ud bs =) “aOaa Sse WY a -d > k- > <t 
Or Cage (SW AD cA NI i lz He AN WIN D Y= - ee 

~ oe Moo Wo. - a ie pe WR ZY OM HR DFO MM WwW to: WY) 
‘so 0CUl J tH ee i - - 2= vs Ww Om eZ NWR OOnyD Ak LJ Oo - 

ji a tum am - em em ttt NEO Zoe) 290. Nr a eo ele 
Ut SS ae isto Te Mt eae! Ze ee JYXINT ZO © me GO © & 
Oo  — Be k- | ot oO ce © SOD w@ a ey xe ANKM Faz ot DS Lee UCU 
re & Trt + aa Ot aA Me MQ Ww HOZO Westie ey et Geter - = 
cc ioe Te Ih > <f sho Wo Zee TF Oe HU ee Qe ela net Saar oS. UW 

t O <4 COOL It <o e— 1D > = oa. 12 oO 9F of = aDWahse —Siae 2. ek OO 
a8 - ba LL or Im ee * AU UUs ob Me UH ORD ck <f WANK OC Zz 
~~ WwW « DOS «w Er = Ss SN YS WL BUR A ZlY IL ID we eK De <Oan~ N & 
re Oe 2 Ja << thie) ea OO f-— MUSH wg iden MA F mia te te Oe 
Ne Ww OO Waey COMIN Wee mer ls w= @ + as DDL ~e NSIS YO > is 
> OO NrEUTt eww 6 wrt ee <tr PUIZZ WZNXitig Ltt tule ew f—--4 OOO € ~ OO 
+ WY & eth rH TS Su Rete OW he OU oO tO Fro ox mime) POM NN Ute tO 
N WwW > «K wrt O felieqT ef WHO MIST D OVOCOONODORIUWYr Ald man LI Om Ol 
~m YF ft aAterr atc KS Om OSORS ADO NZ2AiN earl xreaa SIF Ww OME KH CO HO 
——D Omer Dee Ee eae omy dee) Ee PAZ Ze SIKH NU OOS wR HOO Oe 
as OmMmNayv we te mae @Mnaet WS “~a ew) = FPeHAOMGOrmMme Fe MMYM 
~<T S e ty -~ oom me. om ™ Ny 8 LLILLI ERS uIiaQem LU aw ft] LL t§ OS — ect >< I pe HN C) we (LS 
NOT SBOmOWHROOODRO of oP KINNK CODY ORTON 8 NER ae BW Ne RIUS KOKO 1 mt ee 


SH tH POR Oe OM OK INM HOR OTOWRFOAAE DK OR AU MORT eMOURR Rem LUE MOD SIM 


Zam em eENa eno DOWOWkh-ROTOAHODOUOUNY~OoWwMm << Kem DOINYDW OlVYe OFODVVDe ODRARK INK YK OA AXA 
MOON YR NH bO HOA eH Ste me be et ete tb OY OD BRO me et wer OU wre Ltt eet 
Pa tH ASO eR ee Ok OE IL wm HHO OO Oe Oe MHI IE BI HHO NO & OO Siw feats r ie ey O> O--me 
LIU) meme NNN BPO CS fee Lh LINK pet fe ee pet fe re be ee FR Le RI eh eM hk mb DIO ROA OKO OMT 
MINION ZENA Ow ten oo Ee WM wt Tested Wtwertitides DUI YOODOY HFN NOrRMN~> 
aie Nd Ing | WCOTPOHKIC) RR EIR ER SR ER SE OWOoWOrMi tT OWwWe Sehr O RP DOrRARr=EOD OH HPYOOH + e>HO 
CL Sd ee) OCD wD HE TO Telefe Ra HHA LOY ZIOOK eC HOMY “war Oma ave st NOstesNoonWh a4 
SIN Of Hee ft NS TUAW OR WI CP DYNOCVOYCOCOMOWNY HX INITOUNYOOCOWUC OWA DYOMOOROTKREE HT HK-OOxKS 
KH UIVIVILL Je te NILE de ELL OE NN ote VS BU RUB eb Mt eZ ZOUW CHS ULYHSU CK ELAYNOSULOZOOOUSaS Py woOoOnd 


red EOI AY 


oH ; = HOY 
ms ow onmst W So Oo AN eo N 
© oO mm Om FY v) © Oo ry Mm 
_ 0 OO OY DO Oo 4 wo co Oo oO 


OO OW WOW UWOUOOOW WO WO Wo O UO oO 


46 





© Il 

LJ : 

b~ QO 

= > 

a , 

oO. a 

= ° 

a) wy 

= al 

lk - 

= tit 

oo a ~~ 

VY . i 

ex, lI a 

Lt ' 

ff 3 - 

t 

= Oh - 

“> 

Wi ae ( 

<I >= i 

~ & 
< On OA 
re Hts Wks 
hm ok & 

eS >A >Re 

<< et om 
etx > © ely 
a - em 8 
> WJ m 44 merdinj 
ws oh a be LY ere 
kt >mM a7 
= eC Al ~ ta 
C) mnt cae mH>>< 

= We Lime ry 
IO ~ >< sw 
zwei us (Nt OD we hw 
mt (I QO ~ © @ ‘Ve il nyo 
Lib a oe ot ID Ore oo 

{OO Fl UW ee © wr 
mrad dl HOw DO— Ket Owe 


J. eRe ce 24a 


OOZ~~LOO+~OO0KCON ON cx 


Lh, xo mend Dm Se tee OL SU le TILA 


UOUO 


20 ‘ 
? 


_ 


fi <I & 
NI . 2B > 
: (3 e eam — tml VY) 
© Us = =) < me () 
- Y * €& - + Lu 
j ~ O lL N -_ CoG. 
wd WM et wwe ON ceed ol 
< (omer ase Le > ” LL. 
m= i .- @ os = rot XZ 
Ge -- © > —- ae WV) 
we a) re b- > > + 
| a ieee = | oe ik ete 
WY ~ st Oo 3 ACS 
_ oO XK -~ O > ~ ti<t 
| =e N OC -~ OO 0 NAN. ON 
f+ Oe | AN N O en ie Vv) 
4 re es  & (Coo =a | = Q. 
>e ee N ke Wa <( tu 
-O « > ® ee << ea 
= ~oH Oe ty NN + b- 
(i) Qt Od e os ee — 
= ae (Le Noto b> 4 NO 
~ Oe “~ oC — e~ 8 & 6 = Na re 
oO N ~ WwW NANI FE 2 Cie en) Oo 
mr) ~ II ~ _ -uUrat>S> ll athe SO Ov 
| @ Lu ane LL ee Oe <I<f Or 
J rc) 8)  S: ne Pte ee Lt WiGeae ©) 
Ae ee ee 4 am MR ultb + mor © 
be ~ oO. @ 0 SO Fetttumo - WiC b- 
we cr ~ O Ie nm or+anin amb or 
ee we & b>) fy} MOC OD ge mt 4 EGS ~ oe 2) 
UL ox i wm QO mR toe ORO OD Oo 
C) oH =F & © je ew STUN ODO OLS 
xOY MH Oo) HOY FOmNwww et ONT 
re VY) th. uu C2 Wen © aeseaect eT ~ N> Tre @ = 
a a -~ WY S QD wANb-eHf FRE + 2H F&F O 4 
oe) Ow =- a & Q Re wee eT Os; mH Ue Oo 
4 >N Vv) ri 1 oe et f 
ut “NS eell Ge oH NN LI $e + oH Gia=ige oJ ~— © 
ale =a ON Lj i AS eH DOIHCO-ENN 2& ~Ci-fp- @ LV a 
Ow ftom +O ~ et *ORUNN <clwwNA ew OK ON “ & MO UJ 
iL - De the OD WoO ~ SO RNY POR HMA Om SLUR aN atQ 
le OWO *» I ZO ~— th Were it OR OULIL I > Zao + sem Nee om 
- Ore & Nh - O DO h¥XYA eO>wW HARI QW > ORO J e+ 
> NO ws aN = O SAANHODT JRA teaeaRZamk DE wWeDDi lwo i LU pe 
“ OD Il~}- Wut ve Fe a MAS PWN me HOT ew meet OOD 
eer ~UWi<t pee § TYLL We AIL OID TP ert ed me UNE eet dO et HON #8 
WPOORS §OW OR ODO eee I Om ON ene Sak be em mh OOO I NIT 
Seon adr Oc>S SE ee Wd We ON ee See ZZ He rN LN ee 
TL Mb Nea OO mrt ty St CD) LL, he fe ee be oe Te Oe le DOL OO aw Owe LU 
Tet Se SLL ow ge RCTS Sti co (mari aN oe ee Ne 
et NOS UV 
© et Ov UN ey nd MN 
© © gov NO Oo 0 
0 i Oo O 
> (9k OVO 


U7 





GOO 23 


PISKP ) 


er 
GO 
Lo 


® 
WO 


e 
re 
| 
ae ~ med 
— - ee 
©) A. 
a Q Lis I 
Ud +6 —_ 
OG rt — ~~ ae 
= a i x OC 
~_ ce OQ i) — bh 
_ + +t <[ >> _ 
=> ae oH 3i © - 
+ a. oO Lu ~ cO - 
— as Q >) — @ = 
N rl +t s Li. 5) ae — 
ce © Ww e rt o hr <I bh — 
bk oo] ° —j> = = a WJ OQ 
=> ke Ho Lit Ht Y 3 ~~) Le <[ 
+ = ” Oo - ee = N tL cc 
° co ~ TJ N ~ yy © N oe 
~ ow — OQ <{ ~N WY © — 
~\ ~ © Oo ~— or ~~ D - ae 
~ Y ° <a « e om iu C5 
- — < © “wm * Oo am O _ + 
ale N tt Im OOM to fF Oo cao (5 tu - 
mMO—~ wo am acty oO NO 2 ° Ne © gle - 
OO b- —lU marek NOOO FO Ww O ~~ > ~~ 
Ondiia> ~Y awe Dit JO 48 WO Lore Vn = ~_ h fh Pf -. 

NOOK N 3 pf << NO~waey OAN ao D NO Qaelk eMt ~ Ce Sy a) mmm OC 
Om e@ SIU tr FM ORO ODO Do eM wUMINSEH-O J a en oe SOK ond 
here> << OX ." Hees Cnmaer TM aANNS+ > ~ us e e . - Y~ON x 

ease ee — 3 ee OO CUE OND AM ONMOR ONE QO % 8% vw 4 wos 

ON @ NI enti To aN es OM PN OOYE OM ORS rtiNa my am AANA b= — 8 & 
Om AIR star Oo NORM ORM O OO— © TM Fly ewwn G. mas be eet et ee ma keNe 

>See UO POR OTM UR O | OR WIR Og 02 NI a eee rad le a ed OO oN 

Rewer He COON DHOOM LOX QOVSCTWOMONYK Mok LTO 1OOa0 e ab Oeo. of Os cer OG. NJ ANY 
amo IE NOC me NLU J NR die Fe eam rte IE OD OMY NAO CK OR ARAAD ~—O lib do uuu ee INS Of 
saAOHdO>~OWO;NOK—™ ON DDO ON ODMNMHOrFNYE el @© © ete > |} man meee ~ “NUS 

Oe MD Om HO OH OAK COM CHOY NAR we POS LAs te 8 Nilay yYYtyYYY ee ee 
CNR SOG TW SOM ANE Ht OG HE HCO FD eM HE HEH OR KHANH Ne ee Ye YY XY HW Ne Ge UIOXK 
Limi 35 OCI NR A UROL WIA UNOUWNWONWw EOD POLL RNR FE Fe ORR YY YY VY YY eM S~SDADOZNM 

o> Cet ei | YOO *br- eQO Ch} ss Ww eS DH ete OOS I re tM hhh Oh mw OUI Be oe mem’ 
OCR ADD NOGONAIUETMNAO PNIYUITMAN HOW | LY ICA eOK YY eM Y YY YY Ar DDYOU~KW Owe O~ 
STN> | WRN Em Nw FE RL Rw OR LO HEN KOO Pe OHNO I NW KM www wwew~ HO me rw 
LOON~KES>RHOIHWOUNSH WW ols CYOCU When ma WANK LRH ANNMNO YOM HAW MmwTt 


4b Or HW Ode met i RNR ENN CUI WOM HC KE EE OS I ete te the Er OOF SE eb 


—* 


$) met AY OL DP Sr et te Ne pre et te Ee pe et et ee (QR AIN YOR I Se Y NT UC Th OOO 


LOWE Er eESIOINOO RU TT WOW TOSI NSE Er DENK ECONO ea ca.ca LOQOOLOLCONLO 
MOM et > > > S>OAAKOOKH HYOAOAMHY HLA HOKE SSS raarKKLCAELeYNNNNNUMNNNYZOOVCOSLSUL EL OSNOSIL 


ees a = 
“ MAO No Oo 
0 ONO co 6 oO 
0 NO 


O O OO O 


KeKKK) »SPHI (CK yKKK ) 


Cape 

hm OM Se 
SaSa0M OO 
a i] = “Ll &e 


Serac} 


— _ 
mm ON 
oOo oO 


48 


? 


wv 


— Nc = UNVOMmc& 
— ) 
oe QOOrFO0OCO 
ie ~ ae ae 
Y OOOO ao 
~~ OO Ooo} 
< 
~ on 
~~ ee py, OO, oy py FO 
- ~~ me en (TF) a em 
= el MONOMMAMM 
ae ;I~OQOtkud t¢t 
— GC: CoOcCeaGeay 
- Y ODO oOCeoO 
—J ~ Omineaeeace 
om ~_ OeWnom se 
ww STUVLUN OOM © 
ae eet evete 
a eeeeee¢ee 
we EEE EEE 
x tat J 
— ww vt e@ @« © @ @ @ eCg_n 
eAl fea ee ee 
ery IY ND fF © @ @ @© © © ef 
 fGes “a1 SSO e eae) 
oO we © OC 252 ee 
oO os me eNlVO te <di<i<{<{<) 
= my TLt~ eet eeee 
~< ~~ O. J Us aN me 
~As STEAM MOMOMOYAM 
Com) Cimw HOD OGD oc) 
sod een — Zo. lai © 
<~O-4 NOs 1AODOBOODO } 
7 OMMYe4 mt FE ooo oewe® 
nw J UL LIOOMCMIWMD Bea 
ee xX = 0. SUN ODANOIM SO 
aa Pe Se US Of OP PUVLAN OO fh eS 
asi & OM XM UJoc D+ ee © e 8 ee 
NM SAO YY 1) x NMNeeet*eevweee e 
eA Ilo oe eNuWtiy map ee elULWUUINUIL ee OO Ord eDO OO OO OC WwW 
Om SCOY~DDYD BL FR OOO QO DOU et ort raPrafe rt tet ert ert eo 
he wer ZZ DOTY e ee eee COD He NM A tf A He OFZ 
L<Q PLLA cles Orn "CACY PC HONOKNONHOHONOUNONO wn 
mS OR REE OC ID OC ee ee ee ee we OOF OF OF OF OF OF OF OF OF D 
mm OC — jee Oe OL ret re ee ae St Ne et ee et ed et em tet he 


~Y~DQ OAYCCMOOOKRZ TS WibUuLUuee Tere Lee ret Ome Se res ae 

Zils OMmML OUMNH DOO tt tt tet POD O OO OOAOAUOaA WMA WA OA UMW 
+ = 

o OMe 4s NY MOM Ft WW O ™ DBD WO 

oO cm 0 = 


49 





O 


ele 


12. 


LIST OF (REPERENCES 


Levich, V. G., Physicochemical Hydrodynamics, Prentice-Hall, 1962. 


Medwin, H., “Insitu Acoustic Measurements of Bubble Population in 
Coastal Ocean Waters," Journal of Geophysical Research, v. 75, 
Demos =o1 1, Vantany ae Ce 


Blanchard, D. C. and Woodcock, A. H., ''Bubble Formation and 
Modification in the Sea and Its Meteorological Significance," 
Reuss Ve J, De wte-eocem ay 1d57- 


McCartney, B. S. and Bary sb Vick:, “Bcho-Sounding on Probabie 
Gas Bubbles from the Bottom of Saanich Inlet, British Columbia," 
Deep-séea Research, v. 12, p. 285-294, 1969. 


Garrettson, G. A., "Bubble Transport in the Upper Ocean," paper 
to be submitted to Journal of Fluid Mechanics. 


tase, K. M. and Zweifeie Pr. F., Linear Transport Theory, Addisen- 
Nesiey, i967. in 


Wing, G. M., An Introduction to Transport Theory, Wiley, 1962. 
Wyman, J., Jr., and others, "On the Stability of Gas Bubbles in 
Sea Water,'' Journal of Marine Research, v. ll, p. 47-62, July, 


13 OZ 


LeBlond, P. H., "Gas Diffusion from Ascending Gas Bubbles," 
Veummal of Fluid MNechamiess 72 s5-5p.) (11-710. dee. 


nvestand, F. H., Bubbile Distraputicnmimetic Upper Ocean, Nesteum. 
Thesis, Naval Postgraduate School, 1971. 


Ceschino, F. and Kuntzmann, Jr., Numerical Solutions of Initial 
Value Problems, Prentice-Hall, 1966. 


Haderlie, E. C., Professor, Oceanography, conversation with 
author at Naval Postgraduate School on 18 May 1972. 


20 





ENITIAL DISTRIBUTION LIST 


No.. 


Wenense Documentation Cemmers 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Asst. Professor G. A. Garrettson, Code 61Gr 
Department of Physics 

Naval Postgraduate School 

Monterey, California 93940 


Professor H. Medwin, Code 61Md 
Department of Physics 

Naval Postgraduate School 
Monterey, California 93940 


LT Thomas C. Vayda, USN 


331 W. Ellis Avenue 
Inglewood, California 90302 


oe: 


Copies 


Securntv Classtfication 


DOCUMENT CONTROL DATA-R&D 


(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classified) 





ORIGINATING ACTIVITY (Corporate author) 2a. REF ORT SECURITY CLASSIFICATION 


Naval Postgraduate School Unclassified 


3 REPORT TITLE 


Bubble Distributions in a Quiescent Ocean Calculated from the Bubble 
Transport Equation 


[es REPORT Date 7a. TOTAL NO. OF PAGES 76. NOGOF REES 
f{ June 1972 53 12 


8a. CONTRACT OR GRANT NO. G8. ORIGINATOR’S REPORT NUMBER(S) 


6. PROJECT NO. | 


}92. OTHER REPORT NO(S) (Any other numbers that may be assigned 
thie report) 


‘Tio. O1STRIBUTION STATEMENT 


Approved for public release; distribution unlimited. 


. SUPPLEMENTARY NOTES - . |22. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


13. ABSTRACT 


A bubble transport equation was developed and a formal solution obtained 

| for the bubble distribution function. The distribution function obtained depends 

upon the model used for bubble drag, gas diffusion, and ocean circulation. As a 
specific application, solutions were obtained for a one-dimensional, quiescent, 
steady state ocean with a bottom source at 197 meters. The characteristic 
equations were integrated numerically to test the accuracy of approximate 
analytical solutions. The calculated bubble densities were compared with 
experimental data available for this problem. When gas diffusion is neglected, 


the theoretical and experimental distributions are in general agreement. 


Y 





iS ae ee fe a4 E 5 ST La re cere = we  gaRts 


eo Aa TS 
S/N 0101-807-6811 52 Secunty Classification 25 1ane 








e - Security Classification 


LINK A LINK: B LINK ¢ 


aGiLe 


P@eubbles 

| Bubble Distribution 
Transport Equation | 
Quiesent Ocean 
Bubble Transport 

: 


| 
| ~ 
FO 
dD 1 Ree cal A473 (BACK ) 53 
\ 0 ‘ 5 —) 
101-807-687} Security Classification A-31309 


a 
KEY WORDS 





+m ere “ee 


ed 


ee ee ee 

















wi 


DUDLEY KNOX LIBRARY 








